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CONVERSION  FACTORS,  U.S.  CUSTOMARY  TO  METRIC  (SI)  UNITS  OF  MEASUREMENT 


U.S.  customary  units 
metric  (SI)  units  as 

of  measurement 
follows : 

used 

in  this  report  can  be  converted  to 

Multiply 

by 

To  obtain 

inches 

25.4 

millimeters 

2.54 

centimeters 

square  inches 

6.452 

square  centimeters 

cubic  inches 

16.39 

cubic  centimeters 

feet 

30.48 

centimeters 

0.3048 

meters 

square  feet 

0.0929 

square  meters 

cubic  feet 

0.0283 

cubic  meters 

yards 

0.9144 

meters 

square  yards 

0.836 

square  meters 

cubic  yards 

0.7646 

cubic  meters 

miles 

1.6093 

kilometers 

square  miles 

259.0 

hectares 

knots 

1.852 

kilometers  per  hour 

acres 

0.4047 

hectares 

foot-pounds 

1.3558 

newton  meters 

millibars 

1.0197  x 

10-3 

kilograms  per  square  centimeter 

ounces 

28.35 

grams 

pounds 

453.6 

grams 

0.4536 

kilograms 

ton,  long 

1.0160 

metric  tons 

ton,  short 

0.9072 

metric  tons 

degrees  (angel) 

0.01745 

radians 

Fahrenheit  degrees 

5/9 

Celsius  degrees  or  Kelvins 

To  obtain  Celsius  (C)  temperature  readings  from  Fahrenheit  (F)  readings, 
use  formula:  C  -  (5/9)  (F  -32). 

To  obtain  Kelvin  (K)  readings,  use  formula:  K  ■  (5/9)  (F  -32)  +  273.15. 


5 


SYMBOLS  AND  DEFINITIONS 


D  median  satil  grain  diameter 

d  water  depth 

d^  water  depth  at  seaward  bound  of  shoal  zone 

dg,  water  depth  at  landward  bound  of  shoal  zone 

g  acceleration  of  gravity 

H  wave  height 

Hs  annual  mean  significant  wave  height 

Hs0  ^7  extreme  significant  wave  height  exceeded  12  hours  per  year 
Hs50  annual  median  significant  wave  height 

L  linear  wavelength 

L0  linear  wavelength  in  deep  water 

T  wave  period 

Ts  annual  mean  significant  wave  period 

umax(-d)  maximum  linear  wave-induced  horizontal  velocity  near  the  bed 
x  argument  of  hyperbolic  functions  in  Figure  4 

Y '  ratio  of  density  difference  between  sand  and  fluid  to  fluid  density 

o  annual  standard  deviation  of  significant  wave  height 
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SEAWARD  LIMIT  OF  SIGNIFICANT  SAND  TRANSPORT  BY  WAVES: 

AN  ANNUAL  ZONATION  FOR  SEASONAL  PROFILES 

by 

Robert  J.  Hallermeier 

I .  INTRODUCTION 

Certain  coastal  activities  require  consideration  of  the  seaward  limit  to 
significant  wave-induced  bed  activity  along  a  sand  beach  profile.  Examples  in¬ 
clude  design  of  nearshore  structures,  subaqueous  beach  fill,  and  borrow  or 
disposal  of  material.  The  seaward  limit  to  significant  sand  movement  by  waves 
clearly  depends  on  both  wave  and  sand  characteristics  through  the  mechanisms  of 
wave  agitation  of  sand.  However,  available  guidelines  on  the  seaward  limit 
(U.S.  Army,  Corps  of  Engineers,  Coastal  Engineering  Research  Center,  1977,  pp. 
4-63  to  4-70)  proceed  from  indirect  or  qualitative  consideration  of  wave-sand 
interactions.  This  report  provides  a  new  quantitative  procedure  for  seaward- 
limit  estimation  based  on  definite  wave-sand  interactions.  A  comprehensive 
discussion  of  the  present  procedure  and  of  other  seaward-limit  treatments  is 
provided  in  Hallermeier  (1981). 

The  present  procedure  defines  an  annual  beach  profile  zonation  related  to 
the  seaward  limit.  Figure  1  is  a  schematic  illustration  of  this  zonation  and 
identifies  the  terminology  used  here.  The  water  depths  d^  and  d^,  respec¬ 
tively  the  landward  and  seaward  bounds  of  the  shoal  zone,  are  calculated  in 
Section  II,  using  sand  characteristics  and  statistics  of  annual  wave  climate 
for  a  given  locality.  This  shoal  zone  is  a  buffer  area  where  expected  waves 
have  neither  strong  nor  negligible  effects  on  the  sand  bed  during  a  typical 
annual  cycle  of  wave  action,  according  to  available  knowledge  of  sand  movement 
by  waves. 

The  water  depth,  d^ ,  gives  a  seaward  limit  to  extreme  surf-related  effects, 
so  that  significant  alongshore  transport  and  intense  onshore-offshore  transport 
are  restricted  to  water  depths  less  than  d^.  The  water  depth,  d^  gives  a 
seaward  limit  to  sand  motion  by  usual  waves,  so  that  significant  onshore-offshore 
transport  is  restricted  to  water  depths  less  than  d^.  The  seaward  limit  in 
sandy  regions  reasonably  falls  within  this  calculated  shoal  zone,  although  place¬ 
ment  of  the  limit  depends  on  the  particular  engineering  application,  as  pointed 
out  in  Section  III. 


Figure  1.  Definition  sketch  of  annual  sand  beach  zonation  and  terminology. 


Since  the  shoal  zone  has  a  definite  physical  basis,  the  present  calculation 
procedure  may  reduce  the  need  for  detailed  site-specific  investigations  of  phys¬ 
ical  processes,  and  certainly  aids  in  designing  such  investigations  at  localities 
where  wave  climate  is  known.  Nearshore  measurements  of  wave  climate  along  U.S. 
coasts  provide  a  basis  for  the  example  calculations  given  in  Section  IV.  The 
calculated  shoal  zone  extent  is  to  be  considered  in  conjunction  with  other  in¬ 
formation  pertinent  to  the  seaward  limit  of  the  wave-dominated  sand  profile. 


II.  BASIC  EQUATIONS  AND  SOLUTION  TECHNIQUES 


The  water  depths  d^  and  di  are  calculated  from  relationships  giving  the 
threshold  flow  energy  for  two  distinct  wave  interactions  with  a  sand  bed.  The 
major  assumptions  involved  are  the  accuracy  of  linear  wave  theory  for  calcula¬ 
tions  of  near-bottom  maximum  flow  velocity,  and  a  modified  exponential  distri¬ 
bution  of  nearshore  wave  heights. 

The  landward  bound  to  the  shoal  is  determined  from  the  relationship 


umax(-d) 

Y'gd 


0.03 


(1) 


where 


“maxC-d) 

Y’ 


g 

d 


the  maximum  wave-induced  horizontal  fluid  velocity  near  the  bed 

the  ratio  of  the  density  difference  between  sand  -and  fluid  to 
the  fluid  density 

the  acceleration  due  to  gravity 

the  mean  water  depth. 


Equation  (1)  describes  a  definite  threshold  of  erosive  sand-bed  agitation  by 
steady  wave  action,  and  umax(_<j)  will  be  taken  as  that  due  to  an  extreme  wave 
height  exceeded  12  hours  per  year  in  determining  d^ .  The  water  depth  then  de¬ 
fined  by  equation  (1)  agrees  with  field  data  on  the  maximum  water  depth  for 
significant  profile  changes  (elevation  excursions  greater  than  ±0.5  foot  or 
±0.15  meter)  throughout  a  yearly  cycle  of  wave  climate. 


Linear  wave  theory  gives 

UmaX(-d)  =  T  { sinh  (2ird/L)} 

where  H  is  the  wave  height,  T  the  wave  period,  L 
d ;  and 


L 

tanh (2wd/L) 


(2) 


the  local  wavelength  at 


(3) 


where  L0  is  wavelength  in  deep  water. 

The  modified  exponential  distribution  approximating  typical  measured  near¬ 
shore  wave  heights  is  shown  in  Figure  2.  In  terms  of  standard  annual  wave 
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w? 


nMeosured  Cumulotive  Wove  Height  Data,  22 
Apr.  1948  To  30  Apr.  1949,  Long  Branch,  N.J. 
(Hall  and  Herron,  1950) 


Figure  2.  Graphical  definition  of  statistics  of  exponential 
cumulative  wave  height  distribution. 

statistics,  the  extreme  significant  wave  height  exceeded  12  hours  per  year  (0.137 
percent)  is 

Hs0 .137  =  ^S  +  5-6  0  (4) 

where  Hs  is  annual  average  significant  wave  height,  and  o  is  the  annual 
standard  deviation  of  significant  wave  height.  For  simplicity,  the  period  as¬ 
sociated  withthis  extreme  wave  condition  is  takeft  to  be  the  average  significant 
wave  period  Ts,  with  fair  justification.  The  wave  statistics  Fs,  o,  and 
Ts  should  be  defined  by  one  or  more  full  years  of  daily  or  more  frequent  meas¬ 
urements  . 

Using  equations  (2),  (3),  and  (4),  equation  (1)  can  be  rewritten  as 

I  2vdZ  1  sinh2  l2'ndl.  )  tanh l2lrdZ  )  _  *2(HS  +  5.6  a)2  (5) 

\  L  )  \  L  j  \  L  I  "  0.03y '  (gT|/2ir)i! 

and  the  values  of  y',  Hs,  o,  and  Ts  determine  a  single  value  of  (2Tid^/L).  The 
value  of  dj  can  be  accurately  determined  using  the  iterative  root-finding  cal¬ 
culator  program  presented  in  the  Appendix.  Alternatively,  for  y'  =  1.6  (quartz 
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sand  in  saltwater).  Figure  3  provides  a  graphical  solution  of  equation  (5).  The 
dimensionless  value  of  (2ud^/L)  is  converted  to  d^  by  multiplying  by  {(L0/2ir) 
tanh  (2Trd^/L)}  (see  eq.  3);  Figure  4  provides  the  needed  graph  of  (tanh  x)  ver¬ 
sus  x.  Another  alternative  for  y'  =  1.6  is  the  first-order  analytical  approxi¬ 
mation. 


d,,  =  2HS  +  11  a,  (6) 

which  is  convenient  and  usually  within  ±5  percent  of  an  exact  calculation. 

The  seaward  bound  to  the  shoal  zone  is  determined  from  the  relationship 

2 

umax(-d)  0  ,-,N 

V'gD-  -  8’  (7> 

where  D  is  the  median  sand  diameter  characteristic  of  the  shoal  zone  bed,  e.g., 
at  mean  water  depth  of  (1.5  d^).  Equation  (7)  gives  the  threshold  of  sand  motion 
by  wave  action  in  field  conditions,  and  umax(-d)  will  be  taken  as  that  due  to 
the  median  (50  percent)  wave  condition.  For  this  wave  condition,  equation  (7) 
along  with  the  depth  dependence  in  equation  (2)  defines  a  maximum  water  depth 
for  sand  motion,  d^,  which  corresponds  to  the  seaward  limit  of  the  usual  wave- 
constructed  profile. 

The  modified  exponential  distribution  gives  the  median  significant  wave 
height  as 

HsS0  =  Hs  -  0.3  o  (8) 

The  period  associatedwith  this  median  wave  height  is  also  taken  to  be  the 
average  wave  period  Tg.  With  these  relationships,  it  follows  that 

si„h  (^-) .  | 

The  value  of  dj  can  be  accurately  determined  using  the  calculator  program  pre¬ 
sented  in  the  Appendix.  Alternatively,  after  calculating  the  value  of  the  right 
side  of  equation  (9) ,  the  value  of  (2Ttd j/L)  can  be  determined  from  the  graph  of 
x  versus  (sinh  x)  in  Figure  4,  and  converted  to  dj  by  multiplying  by  {(L0/2ir) 
tanh  (2nd^/L)},  using  the  Figure  4  graph  of  (tanh  x)  versus  x.  Another  alter¬ 
native,  for  y'  =  1.6,  is  to  use  the  empirical  approximation 

d"i  =  (Hs  -  0.3  o)  Ts  (g/5,000  D)0-5  (10) 

which  is  convenient  and  usually  within  ±5  percent  of  an  exact  calculation. 

Table  1  summarizes  the  basic  procedure  for  calculating  the  water  depths  at 
the  two  bounds  of  the  shoal  zone,  for  specified  values  of  g,  y' ,  D,  Hs,  a,  and 
Ts«  The  procedure  can  be  executed  using  either  the  calculator  program  provided 
in  the  Appendix  or  the  graphs  provided  in  Figures  3  and  4.  A  convenient  alter¬ 
native  is  to  use  the  simple  approximations  in  equations  (6)  and  (10) ,  which 


(Hg  -  0.3  o)' 


0.5 


(9) 


8y'gDT 


52 
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Table  1.  Calculation  procedure  for  locating  bounds  of 


the  shoal  zone. 


Bound 

Procedure 

Landward 

1 

Required  site  parameters:  y'  and  5a, 
o,  Ts,  and  g  In  consistent  units. 

Find  dimensionless  root  (2rd{./L)  of 
equation  (5),  either  by  Iteration 
using  the  calculator  program  In  the 
Appendix  or  from  Figure  3  after  cal¬ 
culating  the  value  of  the  abscissa 
(2*  (Hs  +  5.6  o)/gT§}. 

Calculate 

Seaward 

i 

Additional  required  site  parameter: 

D  In  consistent  units  at  water 
depth  of  (1.5  dj). 

Find  the  dimensionless  value  of 
(2rdt/L)  specified  by  equation  (9), 
using  either  the  calculator  program 

In  the  Appendix  or  Figure  A  after 
calculating  the  value  of  the  ordinate 

I”2  <fis  -  0-3  “)2\U-5 

\  8y'gDf2  | 

Calculate 

provide  results  usually  consistent  in  accuracy  with  the  critical  values  given  in 
equations  (1)  and  (7)  to  one  significant  digit.  To  allow  for  possible  effects 
of  tides,  it  is  recommended  that  the  calculated  d^  and  dj  be  taken  as  being 
with  respect  to  mean  lower  low  water  (MLLW)  at  a  locality,  giving  conservative 
(farther  seaward)  locations  for  the  shoal  zone  boundaries.  Also,  it  seems  ap¬ 
propriate  to  use  calculated  results  in  U.S.  customary  units  (feet),  rounded  to 
the  nearest  integer  unit. 

III.  APPLICATIONS  OF  da  AND  dA 

The  shoal  zone  is  an  objectively  defined  buffer  area  of  moderate  wave  effects 
on  the  submerged  beach  profile  throughout  the  year.  Where  wave  and  sand  charac¬ 
teristics  are  known,  the  two  calculated  water  depths  are  useful  in  various 
coastal  engineering  activities  requiring  a  seaward-limit  estimate.  The  possible 
applications  in  the  following  list  are  discussed  in  detail  in  Hallermeier  (1981), 
and  further  evaluations  of  shoal  zone  uses  are  being  conducted. 

(a)  Fathometer  surveys  of  customary  vertical  resolution  (±0.5  foot 
or  ±0.15  meter)  may  not  be  expected  to  resolve  usually  small  bottom 
changes  during  the  year  seaward  of  d£.  Nearshore  surveys  extending 
to  d£  are  needed  to  define  seasonal  conditions  at  a  site. 

(b)  Gages  to  record  waves  incident  on  the  littoral  zone  may  be 
placed  just  seaward  of  the  anticipated  value  of  df ,  so  that  the  meas¬ 
ured  waves  will  not  be  modified  differently  through  the  year  due  to 
significant  changes  in  the  bottom  seaward  of  the  gage  site. 
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(c)  Sand  for  beach  fill  must  be  placed  landward  of  dg  to  in¬ 
crease  the  sand  supply  in  the  intensely  active  littoral  zone.  The 
necessary  quantity  of  nourishment  material  is  at  least  that  required 
to  advance  the  average  profile  landward  of  d^  the  desired  distance. 

(d)  In  calculations  of  shoreline  erosion  due  to  sea  level  rise 
over  long  timespans,  the  ocean  boundary  to  the  equilibrium  nearshore 
profile  may  be  set  at  d^. 

(e)  In  calculations  of  sediment  budgets  pertaining  to  timespans 
less  than  a  few  decades,  the  ocean  boundary  to  the  control  volume  may 
be  set  at  d^. 

(f)  At  sites  with  shore-parallel  contours,  subaqueous  borrow  or 
disposal  of  material  should  be  conducted  well  seaward  of  dj,,  so 
that  destructive  effects  on  nearshore  processes  are  minimized.  Sea¬ 
ward  of  dg ,  borrow  or  disposal  might  be  conducted  in  a  region  with 
irregular  contours  if  bottom  elevations  and  thus  shore  exposure  are 
not  changed  significantly. 

(g)  A  rubble-mound  breakwater  should  be  sited  in  water  deeper  than 
d^  for  the  proximate  region,  if  its  primary  purpose  is  to  provide 
wave  shelter  with  minimum  effect  as  a  littoral  barrier. 

(h)  To  balance  economy  and  function,  it  may  be  advisable  to  build 
a  shore-normal  structure  to  the  length  corresponding  to  the  dj,  con¬ 
tour  for  the  nearby  region  if  its  primary  purpose  is  to  control 
littoral  drift. 


IV.  SHOAL  ZONE  EXTENT  ALONG  U.S.  COASTS 

Tables  2  and  3  present  bounding  water  depths  for  the  shoal  zone  at  selected 
U.S.  coastal  localities.  Using  the  calculator  program  provided  in  the  Appendix, 
calculations  for  the  10  sites  in  Table  2  use  annual  summary  statistics  from 
the  wave  climate  data  reported  by  Thompson  (1977).  For  these  sites,  at  least 
1  full  year  of  nearshore  surface-piercing  gage  data  is  available,  with  summary 
statistics  provided  by  objective  digital  record  analysis. 


Table  2.  Calculated  shoal  zone  bounds  for  10  U.S.  sites  using  annual  statistics 


of  nearshore  wave  gage  measurements  (Thompson,  1977);  y*  •  1.6,  g  •  32.2 
feet  per  second  squared  (1  foot  •  0.305  meter,  1  millimeter  ■  0.0033  foot). 


Site 

Location 

S, 

JliL 

0 

(ft) 

TS 

(s) 

D 

(m.n> 

(ft) 

dl 

(ft) 

Atlantic  coast 

Atlantic  City,  N.J. 

39°21 '  N.,  74°25'  W. 

2.95 

1.55 

8.6 

0.11 

21.5 

82.0 

Virginia  Beach,  Va. 

36°51’  N.,  75°58’  W. 

2.40 

1.40 

8.3 

0.11 

19.0 

65.8 

Nags  Head,  N.C. 

35°55'  N. ,  75°36’  W. 

3.20 

1.80 

8.8 

0.11 

24.2 

88.7 

Atlantic  Beach,  N.C. 

34°43'  N.,  7 6°44'  W. 

2.25 

1.30 

7.2 

0.11 

17.2 

51.0 

Vrlghtavllle  Beach,  N.C. 

34°13'  N  ,  77°47 '  W. 

2.60 

1.10 

7.7 

0.11 

16.3 

66 . 7 

Holden  Beach,  N.C. 

33°55’  K..  78°18’  W. 

2.05 

1.00 

7.5 

0.17 

14.4 

43.2 

Uk*  Worth,  Fla. 

26°37 '  a.,  80°02'  W. 

2.10 

1.15 

6.7 

0.21 

15.4 

34.8 

Culf  of  Mexico  coast 

Naplaa ,  Fla. 

26°08 '  H. ,  $1°49’  V. 

1.10 

0.75 

4.6 

0.12 

9.1 

15.8 

Daatln,  Fla. 

30°  2 3 '  N  ,  86°2V  W. 

1.65 

1.05 

5.7 

0.25 

13.: 

20.'. 

Pacific  coast 

Huntington  Beach,  Calif. 

33°39'  ll..U8o00'  V. 

2.90 

1.05 

13.2 

0.11 

17.8 

139.3 
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Tabic  3.  Calculated  shoal  tone  bounds  for  20  U.S.  sites  using  annual  statistic#  of 
LEO  breaker  observations;  v*  •  1,6,  ft,  ■  32.2  feet  per  second  squared. 


0  -  0.00033 

oot  (1 

foot 

-  0.305 

meter) . 

Sit# 

Location 

1 

in 

0 

n 

n 

ItSl 

(ft) 

GDI 

(ft) 

Atlantic  coast 

kasataague,  rtd . 

3B°U 

75°09' 

W. 

1978 

2.18 

1.14 

17.7 

73.3 

Bull  I, land,  S.C. 

U’55 

N. , 

79°35' 

W. 

1977-78 

Ksa 

0.84 

Bra 

13.7 

52.7 

Tybaa  Lighthouse,  Ca. 

3?°01 

N.  , 

80°50' 

w. 

1976-77 

m 

1.85 

BiSI 

24.3 

8b.  7 

Boca  Raton,  Fla. 

2  j°.!2 

N.  , 

80°04' 

w. 

1971 

1.59 

1.04 

13.7 

31.2 

Gulf  of  Mexico  coast 

St.  Andrews  Park,  Fla. 

30°05 

N., 

B5°40’ 

w. 

1969-70 

®Z1 

1.24 

4.(6 

14.7 

29.6 

Crystal  Baach,  Fla. 

3  )°23 

N.  , 

84°27’ 

w. 

1969-70 

ffifi 

1.42 

4.11 

ES 

27.1 

Gilchrist,  Tex. 

29°31 

N.  , 

94°29' 

w. 

1975 

H3 

EQ 

6a?  5 

Bti 

31.4 

Galveston ,  Tex. 

29°  11 

N. , 

94°58* 

w. 

1975 

1.53 

m 

6.71 

12.5 

39.6 

Corpus  Christ  1,  Tex. 

27°45 

N., 

97°10* 

w. 

1974 

2.59 

1.12 

6.66 

17.0 

6f  .9 

Pacific  coast 

San  Clenente,  Calif. 

33a24 

H., 

117°21# 

w. 

1969 

2.85 

1.40 

15.15 

23.4 

169.4 

Bolsa  Chlca,  Calif. 

33°41 

N.  , 

118b02' 

w. 

1969 

2.42 

1.23 

12.10 

20.0 

114.4 

Pt.  Mugu,  Calif. 

3h"07 

N. , 

119°09' 

w. 

1973 

2.84 

1.00 

14.69 

18.6 

172.1 

Plano  Beach,  Calif. 

35‘  09 

N.  , 

120°39' 

w. 

1969 

3.07 

1.52 

12.03 

24.6 

143.0 

San  Sineon,  Calif. 

35°  34 

N.  , 

121°07’ 

u. 

1969 

3.09 

1.22 

11.58 

21.2 

141.7 

Capitol#  Beach,  Calif. 

36c  39 

N. , 

121°S6* 

w. 

1971 

1.50 

0.87 

11.16 

13.9 

s«.o 

Stinson  Beach,  Calif. 

37°54 

N., 

122°38* 

w. 

1968-69 

3.78 

1.32 

12.62 

23.9 

186.4 

Wright's  Beach,  Calif. 

38°24 

N. , 

123°06* 

w. 

1969 

4.94 

2.10 

11.37 

34.0 

150.2 

Shatter  Cove.  Calif. 

40°  02 

N.  , 

124°  04  * 

w. 

1969-70 

2.33 

1.53 

23.2 

102.1 

Prairie  Creek,  Calif. 

41°21' 

N.  , 

124 ’04' 

u. 

1969-70 

3.44 

1.34 

IPlKg 

22.9 

139.6 

Uapqua,  Oreg. 

43°29' 

M.  , 

124°13 ' 

w. 

1978 

3.74 

1.61 

9.38 

25.6 

172.8 

Calculations  for  the  20  sites  in  Table  3  use  annual  summary  statistics  of 
breaker  observations  collected  under  the  LEO  program  conducted  by  the  U.S.  Army 
Coastal  Engineering  Research  Center  (CERC) .  Such  wave  data  are  liable  to  criti¬ 
cism  regarding  precision  and  subjectivity,  and  seem  to  give  larger  annual  mean 
wave  heights  than  gages  located  seaward  of  the  littoral  zone.  Also,  for  the 
Table  3  calculations,  D  has  been  set  equal  to  0.00033  foot  (0.1  millimeter), 
a  typical  value.  However,  the  Table  3  results  show  consistency  with  Table  2 
results  for  common  regions. 

***************  EXAMPLE  PROBLEM  **************** 

GIVEN:  At  Galveston,  Texas  (29°17’  N.,  9A°47*  W.),  annual  wave  statistics  for 
1966  (Thompson,  1977)  are: 

Mean  significant  height,  Hg  =  1.35  feet 

Standard  deviation  of  height,  o  =  0.9  foot 

Mean  wave  period,  Tg  =  5.7  seconds. 

Median  sediment  diameter  at  the  18-foot  water  depth  off  Bolivar  Peninsula  just 
north  of  Galveston  (U.S.  Army  Engineer  District,  Galveston,  1958)  is: 

D  =  0.00033  foot  (0.1  millimeter). 
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FIND:  The  water  depths  d^  and  d^,  giving  the  extent  of  the  shoal  zone  at 
this  site. 

SOLUTION:  From  equation  (4) , 

Hs0 . 137  =  {1*35  +  (5.6) (0.9) }  =  6.4  feet, 
and  from  equation  (3), 

I  32. 2(5. 7)2  \  ,,,  „  , 

L0  =  | j  =  166.5  feet 

so  that 

H-“— —  =  0.038 
Lo 

and  Figure  3  gives 

27rdp 

-=-=■  =  0.74 

L 

Figure  4  gives 

tanh  (0.74)  =  0.63, 

so  that 

d{.  *  ^ - — ^-2-^|tanh  ^ =  (0.74) (26.5)  (0.63)  =  12.4  feet  =  12  feet. 


Note  that  equation  (6)  gives 

d£  =  { 2 [1 . 35  +  11(0.9)]}  =  12.6  feet  *  13  feet. 
From  equation  (8)  , 

Hs50  =  {1.35  -  0.3(0. 9)}  =  1.08  feet, 

so  equation  (9)  becomes 


sinh 
Figure  4  gives 

and 


'  L  /  1 8(1 . 


TT2  (1.08)2 


6) (32 . 2) (0.00033)  (5.7)' 


0.5 


=  1.61 


/  2irdjN 

sinh"1  (1.61)  =  1.26  =  ( ~~ ) 

/2vd±\ 

{-r)  ■  °-85- 


tanh 
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so  that 


di 


(1.26X26. 45)(0.85)  =  28.33  feet 


Note  that  equation  (10)  gives 


28  feet. 


di  =  { (1.08) (5 . 7) [32 . 25/5, 000(0. 00033) ] 0 - 5 }  =  27.2  feet  =  27  feet. 

With  the  same  input  parameters,  the  calculator  program  in  the  Appendix  gives 
the  more  precise  results: 


12.30  feet  =  12  feet 


d^  =  28.32  feet  =  28  feet. 

*************************************** 


These  calculated  results  are  to  be  considered  in  conjunction  with  other 
possible  indicators  of  the  seaward  limit  to  the  wave-dominated  profile.  Everts 
(1978)  reported  two  geometric  limit  depths  for  Palm  Beach  near  the  center  of 
Galveston  Island  (29°12'  N.,  94°57'  W.).  The  limit  depth  to  shore-parallel 
bathymetry  is  47  feet  (14.3  meters),  and  the  depth  at  the  boundary  between  the 
shoreface  and  ramp  profile  sectors  is  59  feet  (18.0  meters);  both  depths  are 
with  respect  to  mean  low  water  (MLW) . 

Sand-size  variation  along  the  profile  can  indicate  a  seaward  limit  to  wave 
effects  (U.S.  Army,  Corps  of  Engineers,  Coastal  Engineering  Research  Center, 
1977).  For  the  Galveston  area,  Fisher,  et  al.  (1972)  charted  the  boundary  be¬ 
tween  "shoreface  sand  and  muddy  sand"  and  "shelf  mud  and  sand  with  shell"  at 
about  28  to  29  feet  below  mean  sea  level  (MSL) .  (MLW  is  about  1  foot  below 
MSL  at  this  site.) 

Watson  and  Behrens  (1976)  presented  six  profiles  surveyed  to  about  14  feet 
below  MSL  in  June  and  December  1974  and  March  1976  near  New  Corpus  Christi  Pass, 
Texas  (28°4l’  N. ,  97°11'  W.).  Five  of  these  profiles  superpose  to  within  +0.5 
foot  beyond  a  water  depth  of  about  10  to  12  feet  below  MLW.  Charted  profile 
geometry  at  this  site  is  similar  to  that  near  Galveston  (Everts,  1978),  and  wave 
climate  may  be  expected  to  be  similar. 

The  preceding  information  supports  the  calculated  values  of  d^  and  d^  for 
this  site.  In  this  example,  the  shoal  zone  seems  well-founded,  while  the  two 
geometric  indicators  derived  by  Everts  (1978)  are  not  supported  by  other  evi¬ 
dence  on  the  limits  to  wave  effects  on  the  submarine  profile. 

V .  SUMMARY 

This  report  is  a  self-contained  presentation  of  the  calculation  procedure 
(Table  1)  for  locating  a  shoal  zone  along  a  sand  beach  profile.  The  water 
depths  bounding  this  shoal  zone  depend  on  sand  characteristics  (y'»  D)  and  sum¬ 
mary  statistics  of  annual  wave  climate  (Hs,  o,  and  Ts) .  The  shoal  zone  is  de¬ 
fined  so  that  expected  surface  waves  have  only  moderate  effects  on  the  sand 
bottom  in  this  zone  throughout  a  typical  year.  The  calculated  shoal  zone  bounds 
have  applications  in  coastal  engineering  activities  requiring  an  estimate  of  the 
seaward  limit  to  the  wave-dominated  nearshore  region  (see  Sec.  III). 
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The  assumptions  involved  in  the  shoal  zone  model  are  limiting,  to  some  ex¬ 
tent,  and  confirmatory  field  evidence  is  definitely  meager,  but  the  present 
calculation  procedure  is  objective,  and  seems  to  be  a  useful  new  tool  for  engi¬ 
neering  applications.  The  calculated  shoal  zone  supplements  other  estimation 
techniques  for  the  seaward  limit  (U.S.  Army,  Corps  of  Engineers,  Coastal 
Engineering  Research  Center,  1977);  of  these  other  techniques,  analysis  of 
sediment  character  variation  along  the  profile  seems  the  most  defendable. 
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APPENDIX 


CALCULATOR  PROGRAM  FOR  THE  BOUNDS  OF  THE  SHOAL  ZONE 

A  calculator  program  executing  the  Table  1  procedure  is  documented  on  the 
following  four  pages.  Written  for  the  Hewlett-Packard  HP-67  Programmable  Pocket 
Calculator,  this  program  uses  metric  or  U.S.  customary  units,  224  program  steps, 
25  storage  registers,  and  RPN  logic.  The  program  could  be  translated  for  use 
with  calculators  having  other  logic  systems  and  similar  features  and  capacities. 

The  program  incorporates  the  basic  relationships  in  equations  (1)  and  (7), 
the  needed  identities  for  hyperbolic  functions  and  for  linear  waves,  and  an 
effective  interative  root-finding  routine.  The  parameters  to  be  specified  for 
the  calculations  are  Hg,  o,  Ts,  D,  y' ,  and  g.  For  quartz  sand,  y'  =  1.6  in 
seawater  and  y'  =  1.65  in  freshwater,  with  slight  variations  with  water  tem¬ 
perature  and,  of  course,  salinity  (Komar,  1976).  For  most  sites,  especially 
at  sea  level  in  temperate  latitudes,  the  value  of  g  may  be  taken  to  be  32.2 
feet  (981  centimeters)  per  second  squared. 


Program  Deseriplion 


Program  THIe  RoonDinC,  WATER  DepTHS  «F  WAve-BASep  SmOAl  IonE 
Ntm  Robert  4  HallcrmEiER  dm*  Jan.  1981 

Address  RESEARCH  Div«S«ON  ,  U.  S  .  ArmY  CoAVTAl  £NC,tNEe<2tN£-,  l2eSrFAR<rt  C'eNWR 
City  K«NC,mAN  &utwpiM£>,  Foft-T  BelvOiR  Slats  V/kLCsINiA  Zip  Coda  ZIOfeD 


Program  DeacrtpUon.  Equations,  VarteMee.  ate.  The  MAiN  PROGRAM  PERFORMS  Am  ITERATIVE 
SolUTi  #N  Of  TWO  FUNCTIONS  DEFINING  DiHENSiONlESS  WATER.  DEPTH  AT  fAtri 
CND  OF  A  SHOAL  lON  £  ALOM G  THE  SUBMERGED  BEACH  P&OFilE  .  AFTER.  WAVE 
AMP  SAND  CHARACTER  AT  A  SITE  ARE  INPUT,  THE  SECANT  METHOD  FINDS 
A  Boor  accurate  to  three  significant  figures  .  Other  programs  : 
convert  the  dimensionless  boots  to  and  ,  water  depths  in 

METERS  AT  THE  LANDWARD  A*S>  SEAWARD  BOUNDARIES  (  RESPECTIVELY  f 

Perform  units  conversions  ;  and  calculate  critical  velocit \f  w£, 
For  sand  motion  initiation.  Other  sym&ols  are; 

=  RATIO  OF  DENSITY  'J"FfeREWC6‘  BETWEEN  SAND  and  Fluid  TO  Fluid  density 
ANNUAL  AVERAGE  r,»x»NlFlCA>vT  WAVE  PERIOD  )  9  *  ACCELERATION  OF  GRAVITY ; 

Hs»  annual  average  significant  wave  height^  ^characteristic  sand  SitE^ 
or  *  annual  standard  jew*,  n*  '  oF  significant  wave  height^  L»  local  wavelength; 
a.*  near-  Bottom  m>*c*£*ntal  Flow  amPuTude  ;  L^C^T dbet- water Naveiength, 
S*  -  (nr  /l)  »  dimension  less  water  depth;  *  »  dummy  index  standing  For  l  or*.. 


The  two 

o  -  ^(e 


BASIC  FUNCTIONS  ARE  '■ 

•u  _  f«\*  *.  *1 r 

'  0.05 


1  (A^facrl*  i  e-l*  \ 
»»  ' 


G-bl  z-D 


i  +  {o.(e'fc  -  e**)/ (*V  o.5  tr  )  } 


And  o  »  I  T  {a.  (e'5i  -  C**' 

The  secant  method  uses  the  iteration  equation: 

'  %?>  *  %{  -  F«p«W)/CF<*i**  F(V'^)  •  ClBle1 

HER  EQUATIONS  AND  IDENTITIES  EMPLOYED  ARC  :  _ 

*  1YTO./T  »(8  D)ff 
t*wK  Ij*  »  (e>  V* 


asL  si 


Other  equations  and  identTties 
-  He  / 1  si-Tik  ;  1 

siwk  ^  *  i:(cs*  -  e**) 


\0.E 


Clbl  c  1 


UCtawk  L  ;  d*  -(£*  t~k 


a)/^e5a  +  e-?*y. 


Operating  LlmKa  and  Warning*  The  INCORPORATED  BASIC  RESULTS  ARE  VALID 
ONLY  FoR  oPeN- COAST  CONDITIONS  AND  SAND-SllE  SEDIMENTS  (D 
Between  o.oofc  and  o.z  centimeters  ).  Linear  wave  theory 
And  a  modified  exponential  Distribution  For  cumulative  wave 
height  ARE  ASSUMED.  The  Roots  AnD  ^  CAN  SE  PhYSICALLT 
MEANINGFUL  only  IF  they  ARE  BETWEEN  0. 1  AND  S.O  ,  APPROXIMATELY. 


User  Instructions 


/  BOUNDING  WATER  depths  OF 

^  2.A •  iK^wi^d 

r£  SA'.  St*w»r6  k. <A,$i 


t« 

O.Ualat'nv 


WAVE-BASfcDA  SHOAL  TONE 
6E.*U<-t 
roit  Fiidtr 

JcA :  JEc  :k^ 


*► 

i/ 


STEP 

INSTRUCTIONS 

INPUT 

DATA/UNITS 

KEYS 

OUTPUT 

DATA/UMTS 

i 

Load  Program  into  calculator  memory. 

1  II 

,z 

Store  site  parameters  : 

1  II 

A.,  NoRMALCiED  IHM ERSCD  SEDIMENT  DENSiTT 

%‘t  - 

ISTOll  0 

►a.  Annual  average  significant  wwe  Period 

T*.  src. 

ISTOll  1 

Annual  AVERAGE  Si^iJiFltAUT  NAVE  HEIGHT 

H«.  O-L 

Isroll  z 

:i 

fTv  Ft 

|STD||  A 

jUaT.  Annual  STANDARD  DEVIATION  «F  HEIGHT 

<r  .  /*a 

1st®  1  1  4 

\i 

or  r  F*fc 

1  STD  |  1  S' 

j 

W.  ACCELERATION  OF  GRAVITY  ■ 

a  cnw/src* 

1st®  1 1  fe 

o. 

X^t/sec1 

|STO||  7  i 

3 

If  Any  INPUT"  rACAMETEfc  IS  IfsJ  IVON  -  METRIC 

0* 

I  II  1 

UNITS  ,  AND/OR.  OUTPUTS  ARC  WANTED  IN  NcN- 

If  II  A. 

HETRid.  UNITS,  ENGAGE  CONVERSION  PROGRAM- 

1  II  1 

4 

Exchange  yiknwxIsecDHVMix  ecsiSter  ccNTevrs, 
Caul  landward  Bound  FUNCTION  CSa)  • 

If  l|P*S| 

G 

1  2.  II  *  | 

6 

Guess  $a»o.w  and  BE6in  it^Raton  , _ 

1  b  II  1 

1  £  1 1  1 

7 

Wait  for  Stationary  DiS?uAX__Op  fjL*°°r- 
Convert  TO  dj_  IN  METERS  and  STDRt . 

1  II  1 

§cA  - 

8 

If  II  t»  1 

4l( 

0 

Tf  DESiRED,  convert  dx  TO  FCCT. 

If  lie  1 

di.ft 

10 

Exchange  primary  /secondary  register  contents, 
St»RE  additional  site  parameter •" 

If  llrtsl 

M 

1  II  1 

tsi.  Sand  $\te  at  water  Derm  of 

Dj  OWV. 

ISTOll  8  1 

a 

Exchange  Primart/secondarilRKiIStw,  ciWTENts 
Engage  program  to  calculate  and  store- 

IF  IIPRSI 

ib 

1  C  II  I 

CRITICAL  VCUOC  |TY  For  sand  jAETlON  INITIATION  . 

1  II  1 

14 

Wait  fOR  STATIONARY  DISPLAY  OP  CRITICAL  VELaJTY 

1  II  1 

NZ,  'Ws« 

i  S 

Call  seaward  bound  Function  ( $*■) . 

IS  II A  1 

lb 

Guess  I.Z  AnD  BEGIN  ITERATION. 

1  1  1 1  1 

|  7.  1  1  1 

|ej|  | 

»7 

Wait  FOR  STATIONARY  display  of  pcot 

1  II  1 

f;  — 

18 

Convert  ^i.to  d;  in  meters  and  stor£. 

IF  lid  1 

dA,  >wv 

•9 

If  desired  ,  convert  d;  To  Farr 

IF  lie  1 

d^  Ft 

zo 

foR  new  CASE,  Clear  All  registers 

IF  IM* 

AND  RETURN  to  STCP  Z, 

IF  )l?*si 

IF  IMa* 

2  1 


Program  Listing 


KEY  ENTRY 
e\  Lbt 


0 

00 

. 

8  S 

♦ 

04 

6 

08 

STOO 

55  O') 

Rcl  S 

54  05 

C*30 

51  it 

C.ToS 

22  OS 

F7o 

51  82 

•STo  L 

55  02 

f  i-Si.  S 

>1  IS  0 S 

iZclO 

55  00 

RCL  5 

54  OS 

_ COMMENTS 

Omits  C0M»er.smM 

AMD  STDRA&C. 


STEP  KEY  ENTRY 

PSP  4 

_  f  CaSft  8 

_ gO.  B 

°«>  faro  o 

_ F  MSI  b 

_ gCL  0 

_  I  6SK2) 

_ STO  ft 

_  F  U5l  0 

_ get  » 

_ RCu  o 

_ STD  A 


KEY  CODE 
25  OS 
St  22  II 

54  12 
22.  00 

Si  2g  Ob 
>4  00 

?'  M. 

55  12 
SI  2.?  00 


SI  SI 
22  07 
SI  82 
SS  04 


St  2S  07 


Stoke  Function 
nuhEEE. 


FiNttC  -  DiFf6*.e«ce 
AVftOKtMATlON  FOR. 
first  oeRiYATivE. 


5S  22 _ 


51  25  1 5  I  Calculate-  and 
STDRf  CRITICAL 
4  00  I  YELDCiTT. 


MFNSlOMuCSS 

lAmpyyaRD  Sound 
poMcnoiu. 


|  *CD  »  • 

TJ  O  « 

V  Q.3  ® 

«  *M  «  OO  3 
I  OC  'J  C 

>  c  •*  «n  ^  o 

►  **4  >  C*  U  C 

4  u  h  ^  •  H 

(  W*  <?>  S 


s 

u  « 

*s 

x  ^  S 
o*u  a. 

*s  « 

•  c  o> 

o 

*  l~ 

«a 

■H  ^*3 

•  ^ 

6  S  c 

SUV 

o  ai  o 

U  H  N  C 

I  «m  w  |  01 

l  -•  u 

1  .  CO  4> 

i  0J 

r«  «oi 

•S  o  l. 

>  *  s  c 

H  B  h 

•2  °  ~  « 

,  s  o 

i  2  r-.  t4  ^4 

!  ®  n  oi  ^ 

>  4J  CL 

..  c  <0 

& 

*5  ^  ^  ° 

*>  H  y  H 

^  ^  L  -O 

|  *QO  ..  0)  -O 

tj  cr>  so 

4  «M  V  dp  9 

i  M  O  — 

>  c  -H  n  ^  u 

>  «H  >  CM  K  C 

4  V4  l*  < — •  V  *-* 


H  U  ••  C 
«  H  O 

3  O  --4 

5  ®  3 

•  i  s  § 

s.a| 

oi  *-< 

«?«-2 

0>  *  W  « 

>  »  «  O 

SH  0)  ^ 

■2*5 

x  a:  a!  o 

-  .  s 

■U  *'>  flO 

l>  C  rH 

ist  a 

S  V  V  o 

ei  «•< 

2  S  5  2 

•o  i‘“  , 

S  * 

S"m  a 

S  3  5  • 

«  «H  <0  ~4 

Z  ■*  o  ji 

X  ^  «  « 

tti  O  H 

“  S.  2 
§  <«  S 
•1  = 
o  S 

J  i 

^  ■ 

« 

L  C  *  ^4 

5  0  K  *4 
<H  *4  00 

m  «  o  c  - 

*  3  52 

8  o8  tn  i 


r-i  *J  *•  c 

to  v  o 

3  0  *r« 

C  *•  u  *J 


*  S  U  -H 

(I  L  « 

>  o»  w  o 

!h!| 

Jas 


^  UH 
•  *4  U  A 

oo  ..  «  ^ 


*2 

V  00  3 
O  — ^  C  «H 
HWf  U 
>  N  h  C 

M  «— i  V  *H 


5  y 

,  >  r*.  n  «h 

,  *J  O. 

..  ..  c  to 

a  & 

i 

»  2  u  2 
-•  <0 
•  00  .«  01-0 
tl  ON  « 

®“  i*S 

•  «4  *  *■  TJ 

i  t*4  oi  op  5 

t  OU«fiH 
>  c  ®4  n  ■»<  u 
i  M  >  N  M  C 

«  hi  w  —  41  h 

I « <X  s 


